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We study the high temperature magnetic order in SrRu 2 06 by measuring magnetization and 
neutron powder diffraction with both polarized and unpolarized neutrons. SrRu 206 crystallizes 
into the hexagonal lead antimonate (PbSb206, space group P31m) structure with layers of edge¬ 
sharing RuOe octahedra separated by Sr^"*" ions. SrRu206 is found to order at rjv=565K with Ru 
moments coupled antiferromagnetically both in-plane and out-of-plane. The magnetic moment is 
1.30(2) ^b/Ru at room temperature and is along the crystallographic c-axis in the G-type magnetic 
structure. We perform density functional calculations with constrained RPA to obtain the electronic 
structure and effective intra- and inter-orbital interaction parameters. The projected density of 
states show strong hybridization between Ru 4ci and O 2p. By downfolding to the target t 2 g bands 
we extract the effective magnetic Hamiltonian and perform Monte Carlo simulations to determine 
the transition temperature as a function of inter- and intra-plane couplings. We find a weak inter¬ 
plane coupling, 3% of the strong intra-plane coupling, permits three-dimensional magnetic order at 
the observed Tjv. 

PACS numbers: 75.50.Ee, 75.47.Lx,75.30.Et,71.15.Mb, 


It is evident that for the design of the next genera¬ 
tion of multifunctional devices, we need new paradigms, 
new principles and new classes of materials. Magnetism 
is arguably the most technologically important property 
arising from electron interactions. One of the central 
questions has been how to create magnetic materials with 
high transition temperatures Tc for room temperature de¬ 
vices. Broadly two paradigms define the formation of the 
magnetic state starting with fermions at finite tempera¬ 
tures: for weak Coulomb interactions U/W compared to 
the bandwidth, one expects a Fermi liquid at finite tem¬ 
peratures, followed by a Fermi surface nesting instability 
that opens a gap in the spectrum resulting in a Slater 
antiferromagnet below Tc ~ W which is expo¬ 
nentially suppressed in the coupling (c is a constant). In 
the opposite regime for \JjW I, local moments form 
on a much higher temperature scale T* Ki U opening a 
large Mott gap and order on the scale of antiferromag¬ 
netic (AF) superexchange J ^ W'^/U. The transition 
temperature as a function oi U/W reaches its maximum 
in the fluctuating regime between these two regimes. 
SrTcOa is currently believed to be at the maximum with 
a Tc ~ 1000 K [T}^. Tuning U/W by combining differ¬ 
ent 3d and 5d ions in double perovskites also increases Tc 
well above room temperature [7], as in Sr2CrRe06 with 
a,Tc = 635 K [S] and Sr2Cr0s06 with a = 720 K [5]. 
These are examples of Mott-Hubbard antiferromagnetic 
insulators. NaOsOa orders at Tc = 410 K nnmi], which 


is a rare example of Slater insulator approaching the fluc¬ 
tuating region from the itinerant side. All the above 
reported high Tc compounds have perovskite structures 
with a d® electronic configuration. 

Recently, a metastable compound SrRu206 with a 
quasi two dimensional structure was proposed to or¬ 
der antiferromagnetically with Tjv above 500 Knu. 
SrRu 2 06 crystallizes into the hexagonal lead antimonate 
(PbSb206) structure with the space group P31m. As 
shown in the inset of Fig. 1, the structure consists of lay¬ 
ers of edge-sharing RuOg octahedra separated by Sr^+ 
ions sitting in the oxygen octahedral interstices. In the 
ab plane, the Ru ions form a honeycomb array. In the 
previous study [T2], the magnetization was measured up 
to 500 K without finding a signature for the magnetic 
transition, though the room temperature neutron powder 
diffraction observed extra reflections absent in x-ray mea¬ 
surements. In this Letter, we report our magnetic and 
neutron diffraction study of SrRu206 up to 750 K. Our 
diffraction measurements with both polarized and unpo¬ 
larized neutrons confirm that SrRu 206 orders antiferro¬ 
magnetically at Tn = 565 K with a magnetic moment of 
1.30(2) ^b/Ru along c-axis. The magnetic measurements 
suggest that strong two-dimensional magnetic correla¬ 
tions persist above Tjy as highlighted in Fig. 1 above the 
kink. Our band structure calculations show strong Ru 4d 
and O 2p hybridization, which accounts for the reduced 
moment of Ru®®" ions and possibly the high magnetic or- 
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FIG. 1. (color online) Magnetic susceptibility in the temper¬ 
ature range 300 K < T < 750 K measured upon cooling in an 
applied field of 50kOe. The solid curves highlight the slope 
change at 565 K. Inset shows the crystal structure. 


dering temperature as in SrTcOa. We analyze the differ¬ 
ent exchange pathways and derive an effective magnetic 
model with predominantly antiferromagnetic interactions 
within the plane, with exchange parameters determined 
from downfolding and constrained RPA, as well as anti¬ 
ferromagnetic interactions between planes. Using Monte 
Carlo simulations of such a coupled layer model we cal¬ 
culate the dependence of T^r on interlayer coupling and 
by comparison with the experimental T/y determine the 
interlayer magnetic interaction to be small, about 3% of 
the in-plane interaction. 

We synthesize polycrystalline SrRu206 by a hydrother¬ 
mal technique as reported previously [T2] and confirm 
that the obtained powder is a single phase from room 
temperature x-ray powder diffraction. All the peaks ob¬ 
served can be indexed with the space group P31m and 
lattice parameters 0 = 5.2037(2) A and c = 5.2341(4) A. 
Fig. 1 shows the temperature dependence of magnetic 
susceptibility measured in the temperature range 300 K 
< T < 750 K using a Quantum Design magnetic property 
measurement system. The data collected in warming and 
cooling processes overlap suggesting little or no sample 
decomposition below 750 K. Above room temperature, 
the magnetic susceptibility increases with increasing tem¬ 
perature. As highlighted by the solid curves, there is a 
slope change at 565 K, which signals possible long range 
magnetic order. 

To study the nature of this slope change, we carry out 
neutron diffraction experiments on about 0.4 g powder 
sample in the temperature range 40 K < T < 600 K using 
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FIG. 2. (color online) (a) Room temperature x-ray and neu¬ 
tron powder diffraction patterns in a narrow d-range high¬ 
lighting the extra reflections observed by neutron diffraction, 
(b) and (c) show the (1 0 0.5) magnetic and (1 0 1) nuclear 
peaks measured with the polarized neutron in the horizontal 
field configuration Pq || Q at room temperature, respectively. 
The observed weak intensity in (c) with flipper on comes from 
the finite instrumental flipping ratio which we estimate to be 
1/10 by comparing the integrated intensity of the (—h) and 
(-I--I-) scans of (1 0 1). (d) The temperature dependence of 
integrated intensity of (1 0 0.5) and (1 0 1) peaks. The solid 
curves are a guide to the eye. Inset shows the G-type mag¬ 
netic structure. 


the HB-IA and HBl triple axis spectrometers and HB- 
3A four-circle spectrometer located at the High Flux Iso¬ 
tope Reactor at Oak Ridge National Laboratory [T3]. As 
shown in Fig. 2(a), neutron powder diffraction observed 
reflections that are absent in x-ray powder diffraction 
pattern. To confirm the magnetic origin of these extra 
reflections, we performed diffraction measurement using 
polarized neutrons with the neutron polarization parallel 
to the moment transfer: Pq || Q. With the spin flipper off 
or on, we measured both the (H—h) non-spin-flip and the 
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(—h) spin-flip scattering of different reflections of inter¬ 
est. Fig. 2(b) and (c) show the intensity of (1 0 0.5) and (1 
0 1), respectively. As discussed in Refs. [I4HI6] . coherent 
nuclear scattering is always (-I--I-) non-spin-flip scattering 
because it never causes a reversal or spin flip of the neu¬ 
tron spin direction upon scattering. On the other hand, 
magnetic scattering depends on the relative orientation 
of the neutron polarization Pq and the scattering vector 
Q. Only those spin components which are perpendicular 
to the scattering vector are effective. Thus for a fully 
polarized neutron beam with the horizontal held config¬ 
uration, Pq II Q, all magnetic scattering is (—h) spin-flip 
scattering, and ideally no (-1—1-) non-spin-flip scattering 
will be observed. The strong intensity shown in Fig. 2(c) 
was observed in the (-I--I-) channel as expected for the nu¬ 
clear (1 0 1) peak. The strong scattering detected in the 
(—h) spin-flip channel confirmed the magnetic origin of 
the (1 0 0.5) peak. 

The diffraction study with polarized neutrons clearly 
shows that the extra reflections come from a long range 
magnetic order instead of any structural transition. 
Those extra reflections can be indexed on the basis of the 
magnetic scattering with a G-type AF structure with the 
magnetic moment direction along the crystallographic c- 
axis. The Rietveld refinement of the neutron diffraction 
patterns m yields a magnetic moment of 1.30(2) /tb/Ru 
at room temperature. The magnetic moment increases 
slightly to 1.34(3) /tb/Ru upon cooling to 40 K. The mo¬ 
ment is smaller than the spin moment oi 3 fj-B expected 
for a a half-filled t 2 g band, which signals a strong cova¬ 
lency of the Ru-0 bonds. 

We further follow the temperature dependence of the 
integrated intensity of the (1 0 0.5) magnetic peak in the 
temperature range of 40K < T < 600K. As shown in 
Fig. 2(d), the (1 0 0.5) magnetic peak disappears around 
565 K, where a slope change is observed in Fig. 1 in the 
temperature dependence of the magnetic susceptibility. 
Figure 2(d) also shows the evolution with temperature of 
the integrated intensity of the (10 1) nuclear peak, which 
shows little temperature dependence in the temperature 
ranged studied. The coincidence of the disappearance of 
(1 0 0.5) magnetic peak and the slope change in magnetic 
susceptibility suggests that a G-type long range magnetic 
order takes place at 565 K in SrRu206. 

In SrRu206 the Ru ions are in a (P electronic configu¬ 
ration in edge-shared octahedral cages formed by the O 
atoms. It is important to note that the Ru-O-Ru bond 
angle is close to 90° hence both ferromagnetic (F) and 
AF mechanisms through intermediate oxygens are ac¬ 
tive in fourth order processes according to Goodenough- 
Kanamori rules [niiis]. There are three competing pro¬ 
cesses that contribute to the exchange interaction: (a) 
the direct overlap of the half filled t 2 g orbitals produces 
a second order AF interaction, (b) The transfer of elec¬ 
trons between an oxygen pz orbital and Ru dzx and dyz 
orbitals on two neighboring Ru atoms results in an AF 


superexchange coupling, (c) The transfer of electrons be¬ 
tween Ru t 2 g orbitals and mutually orthogonal oxygen p 
orbitals results in a F interaction driven by Hund’s cou¬ 
pling on oxygen. 

It is rather intriguing that SrRu206 orders at such a 
high temperature, given the competing magnetic inter¬ 
actions. To estimate the relative magnitude of the above 
competing interactions and to understand the mecha¬ 
nism for the high AF ordering temperature, we per¬ 
form DFT calculations including constrained RPA to ob¬ 
tain effective hopping and interaction parameters. We 
also perform spin density functional calculation with 
the WIEn2k [19] package using the exchange-correlation 
functional proposed by Perdew et al. [20] and the full- 
potential linearized augmented plane-wave method in¬ 
cluding the spin-orbit coupling. In the calculation, we 
use the experimental lattice parameters and atomic con¬ 
figurations determined at room temperature [12]. 

Our DFT calculation shows that the G-type AF state 
is the most stable. The non-magnetic (NM) and the C- 
type AF (i.e., AF in the a6-plane and F along the c-axis) 
states have higher energy than that for the G-type AF 
state, and there is no ferromagnetic or A-type AF (i.e., 
F in the a6-plane and AF along the c-axis) metastable 
solution. These results indicate that there is a strong in¬ 
plane AF correlation compared with that for the out-of- 
plane direction. Galculated local spin magnetic moment 
of Ru atoms is about 0.9 /iB per atom for the G-type AF 
state and in reasonable agreement with our experiment. 
Fig. [^shows the band dispersion and (projected) density 
of states (DOS) for the NM state. The Og bands [23] 
are well isolated from the Cg and oxygen p bands. Strong 
Ru-0 hybridization is seen in the DOS, which points to 
the origin of the reduced local spin magnetic moment of 
Ru atoms. 



FIG. 3. (color online): (a) Band dispersion and (b) (pro¬ 
jected) density of states for the non-magnetic state. 

Using a density response code m recently developed 
for the Elk branch of the original exciting fp-lapw 
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code im, we derive low-energy effective models. Start¬ 
ing with the DFT calculation for the NM state, we first 
construct the Wannier functions and calculate the trans¬ 
fer integrals between them. Next, we evaluate the inter¬ 
action parameters, the Coulomb repulsion U{r) and the 
exchange coupling J (r) by the constrained RPA . For 
simplicity, we neglect the spin-orbit interaction since it 
is irrelevant in the following analysis. 

To evaluate the AF coupling through the direct over¬ 
lap between d orbitals, we derive a low-energy effective 
model for the Ru t 2 g and O p bands in the energy window 
[—7.0:-|-1.0] eV. We obtain the onsite Hubbard Ud = 5.3 
eV and the largest transfer hopping between d orbitals 
tdd = 0.19 eV, which result in a small AF coupling, 
J - ^tlJUd = 0.03 eV. 

We also derive an effective model only for the Ru 
t 2 g bands in the energy window [—1.4:-|-1.0] eV. In this 
model, the Wannier functions are the Ru t 2 g orbitals hy¬ 
bridized with the surrounding O p orbitals, which allows 
a direct evaluation of the superexchange couplings. We 
obtain the on-site Hubbard U = U{r — 0) = 2.7 eV, the 
Hund’s coupling — J{r = Q) — 0.28 eV, the nearest- 
neighbor off-site Coulomb interaction V = 1.1 eV, and 
the largest nearest-neighbor transfer hopping t = 0.28 
eV. For the AF superexchange coupling, these values re¬ 
sult in Jaf 4t^/([/ — V) = 0.20 eV. On the other 
hand, the F superexchange coupling Jp is evaluated as 
the nearest-neighbor off-site direct exchange ~ 0.03 eV. 
The superexchange AF coupling Jaf dominates over J 
and Jf- These estimates are replaced by an exact treat¬ 
ment in the following analysis. 

In the atomic limit, each site has an S = 3/2 spin. 
When the two sites are coupled, the eigenstates can be 
labeled by the total spin ^ = 0,1, 2, 3. Using the effective 
values of t, U, and Jh in the t 2 g effective model [13], 
we perform exact diagonalization for two sites to obtain 
the energies of states labeled by S. These eigenvalues 
determine the exchange constants of a general effective 
spin-3/2 Hamiltonian for two sites as Heff = Eo + Ji{Si ■ 
S 2 ) + J 2 {Si ■ 82 ^ + JsiSi ■ 82 )^. We find Ji = 45.6 meV, 
J 2 = —2.0 meV, and J 3 = 0.5 meV. The interaction at 
the two site level is then primarily Heisenberg AF. 

For purposes of modeling the system, we retain the 
Heisenberg AF nearest neighbor interactions within the 
plane (Jy, 3 neighbors) taken as Ji. Also, motivated 
by the experimental observation of G-type ordering, we 
introduce, in addition, a coupling between nearest neigh¬ 
bor planes (Jl, 2 neighbors). The classical Hamiltonian 
describing magnetism in SrRu 2 06 is given by: 

T^classical — ^ ^ 8i ' 8j J ^ ^ 8i ■ 8j (1) 

(bill (b>-L 

For J± = 0, according to the Mermin-Wagner theorem 
the long wave length spin waves destroy the long range 
order and consequently, the transition temperature tends 
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FIG. 4. (color online). Classical Monte Carlo simulations 
using the model of Q showing the staggered magnetization 
per site versus temperature in units of Jy for (a) Jx/Jy = 0.03 
and (b) J±/J\\ = 10“^ for three system sizes L (total system 
size 4L®). Lowering J± at fixed Jy reduces the transition 
temperature. The transition temperatures extracted by finite 
size scaling are O.lOJy and O.SOJy respectively. 

to zero in the thermodynamic limit. However, a small 
inter-plane coupling can stabilize magnetic order. Our 
DFT calculations show that the inter-plane hopping pa¬ 
rameters are small with respect to the in-plane parame¬ 
ters (tx/t|| ~ 0.1) which suggest that the AF inter-plane 
couplings J± are between 0.1 Jy and O.OlJy. We perform 
classical Monte Carlo simulations on a layered honey¬ 
comb lattice to obtain the transition temperature as a 
function of the ratio of these coupling constants. Fig.|^ 
gives the results of these simulations for two values of 
Jx/Jy. 

The transition temperature obtained from Monte 
Carlo Tj^®(Jy,Jx) depends on the in-plane and out- 
of-plane couplings in general as seen in Fig.|^ Using 
our estimated value for Ji obtained from downfolding 
and exact diagonalization, we obtain the AF Heisenberg 
spin-3/2 coupling constant Jy = (3/2)^ Ji = 102.6 meV 
or 1190 K. Further by using the experimental transition 
temperature of 565 K, we find that Jx = 36 K or equiv¬ 
alently J±/J\\ « 0.03 fits the experimental results. The 
theoretical results also suggest that the transition tem¬ 
perature can be enhanced by increasing the in-plane cou¬ 
pling, for example by chemical or applied pressure. 

As illustrated in the inset of Fig. 1, the quasi two di¬ 
mensional crystal structure of SrRu 206 distinguishes it¬ 
self from other reported high Tc compounds with a per- 
ovskite structure. With the nonmagnetic Sr spacing lay¬ 
ers the inter-plane coupling is expected to be weak. This 
is supported by our Monte Carlo simulations. The inter¬ 
plane coupling is small, but critical for the three dimen¬ 
sional magnetic order. The strong in-plane magnetic in¬ 
teraction is dominated by the AF superexchange cou¬ 
pling between rutheniums mediated by oxygen. We no¬ 
tice in Fig. 1 only one weak slope change around Tjv and 
also the magnetic susceptibility increases linearly with in¬ 
creasing temperature above T^. The absence of a Cure- 
Weiss-like paramagnetic behavior above suggests that 
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strong two dimensional magnetic fluctuations exist above 
Tat and persist until this compound decomposes around 
800 K. 

Our theoretical modeling highlights the mechanism for 
the large AF temperature. Similar to the mechanism for 
the perovskite SrTcOa in Ref. [2], we find that SrRu 2 06 
is close to the fluctuating regime with U^W on the local¬ 
ized side. This is facilitated by substantial hybridization 
between Ru 4d and O 2p within the unit cell as seen 
from the projected density of states in Fig. 3(b) that re¬ 
duces the effective U. In addition there is large hybridiza¬ 
tion between unit cells generating a large band width VF, 
bringing this material close to the crossover region where 
Tat is enhanced. Thus, the large covalency reduces the 
Ru moment but also facilitates a strong (in-plane) AF 
superexchange interaction important for the high Tpf. 

With multiple t 2 g orbitals and Coulomb correlations 
that generate magnetism on a honeycomb lattice, it is 
possible that this material has interesting topological 
properties that still need to be explored, especially as 
we replace Ru with the heavier Os and spin-orbit cou¬ 
pling becomes important. SrRu 2 06, therefore, provides 
a new materials platform for studying the mechanism in¬ 
ducing high temperature magnetic order and other exotic 
phenomena in Ad and 5d transition metal oxides. 
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